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FOREWORD 

This report was prepared by Lockheed's Huntsville 
Research & Engineering Center, under Contract NAS8-26338, 
"Response of Aerospace Structures to Reentry Aerodynamic 
Excitation, " for the George C. Marshall Space Flight Center 
of the National Aeronautics and Space Administration. The 
period of performance was from 26 June 1970 to 25 October 
1971. The work was administered under the technical di- 
rection of the Astronautics Laboratory with J.B. Herring 
acting as the Contracting Officer's Representative. 
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SUMMARY 

An extensive literature survey was conducted of aerodynamic excitation 
and associated structural response. A bibliography based on the survey is 
presented. Very little applicable data under actual reentry conditions was 
found. The survey did reveal many wind tunnel tests and some in-flight data. 

A study was made of aerodynamic excitation and proposed space shuttle 
configurations and reentry trajectories. The study revealed that attached 
turbulent boundary layer, separated turbulent boundary layer and base pres- 
sure fluctuations are the regions of aerodynamic excitation most applicable 
to exterior panels of space shuttle -type vehicles. Prediction techniques are 
presented for these three aerodynamic environments. 

An investigation was made of methods for predicting the response of 
panel structures to acoustic excitation. A normal mode approach was developed 
for the low frequency response and a statistical energy approach was utilized 
for the high frequency response. 

A finite element model was made of a thermal protective system (TPS) 
test panel. The Structural Network Analysis Program (SNAP) was used to 
determine the natural frequencies and mode shapes of the panel. A computer 
program was generated to calculate the response of the panel to a reverberant 
acoustic field, using the normal mode approach developed for low frequency 
response. The power spectral densities of the response displacement and 
acceleration are presented. 
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Section 1 
INTRODUCTION 

With the advent of the space shuttle concept for manned space flight, 
a new regime of aerodynamic excitation has evolved, that of the reentry noise 
environment of a lifting body. In order to cope with this new environment, 
information was needed on the nature and magnitude of this excitation and 
the expected response of space shuttle-type structure to this excitation. To 
help supply this information, Lockheed undertook this study of the response of 
aerospace structures to reentry aerodynamic excitation. 

An extensive literature survey was made of aerodynamic noise data and 
associated structural response. A listing of the applicable documents reviewed 
is given in the Bibliography. The survey revealed very little usable data under 
actual reentry conditions. Consequently, the prediction techniques developed 
for reentry aerodynamic excitation are based primarily on the results of wind 
tunnel tests and flight measurements taken on airplanes. Analytical methods 
were then developed to calculate the expected response of the external panels 
of space shuttle -type vehicles to this excitation. The methods developed are 
also applicable for predicting the response of the structure to the expected 
acoustic excitation during liftoff and the response of test panels of the proposed 
exterior surface to reverberation chamber tests. The response analysis 
methods were used to predict the response of a Haynes 25 alloy L-605 thermal 
protection system test panel in a reverberation chamber test. 
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Section 2 

REENTRY AERODYNAMIC EXCITATION 

An extensive literature review of aerodynamic noise data and analyses of 
proposed space shuttle configurations and reentry trajectories revealed that the 
primary regions of aerodynamic excitation applicable to panels on the space 
shuttle are attached turbulent boundary layer, separated turbulent boundary 
layer and base pressure fluctuations. The attached turbulent boundary layer 
excitation is applicable to the bottom surfaces of the shuttle at the high angles 
of attack associated with high Mach numbers and to most protube ranee -free 
surfaces at low angles of attack associated with low Mach numbers. Panels 
immediately in front of contour discontinuities and protuberances are excited 
by separated turbulent boundary layer. The base pressure fluctuations are 
applicable to the top surface of the shuttle at the high angles of attack associ- 
ated with high Mach numbers. 

Of these three excitation sources, the attached turbulent boundary layer 
has been the subject of the most extensive studies. The studies concerning sep- 
arated turbulent boundary layer and base pressure fluctuations are less numer- 
ous and less detailed. In all cases, the primary sources of reliable data are 
wind tunnel tests. Most flight data suffer from instrumentation problems and/or 
unknown flow conditions and parameters. Consequently, in developing the pre- 
diction technique for the various flow conditions, flight measurements were 
used only in determining the overall fluctuating pressure levels. The power 
spectra, correlation functions and convection velocities were obtained from 
wind tunnel tests. 

2.1 ATTACHED TURBULENT BOUNDARY LAYER 

An analysis was made of the attached turbulent boundary layer. A de- 
tailed description of the boundary layer and the problems encountered in measur- 
ing boundary layer parameters and fluctuating pressures is given in Appendix A. 
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The literature review revealed considerable scatter in the measured data of 
attached turbulent boundary layers. Possible sources of this scatter include 
differences in microphone size, microphone mounting and surface roughness, 
wind tunnel noise and variations in flow parameters. Due to these wide varia- 
tions in the data and the lack of data in certain Mach number ranges, a purely 
empirical approach could not be used for predicting the bound ary -layer excita- 
tion. A combination of theoretical and empirical techniques was used. 

The prediction techniques for the attached turbulent boundary layer were 
taken from the studies of Houbolt (Ref. 1), Blake (Ref. 2), Bull (Ref. 3) and Bies 
(Ref. 4). Houbolt was selected for this theory on the effect of Mach number on 
the overall fluctuating pressure level. Blake was selected because of his superior 
test conditions and instrumentation. Blake's experiments were conducted in the 
low turbulence acoustic wind tunnel at MIT and he used pinhole microphones to 
provide high frequency data without extensive transducer- size corrections. In 
many areas, Blake's results were in close agreement with previous investiga- 
tors. Bull was selected as one of the earlier investigations whose results were 
verified by Blake and others. Bies made an empirical study of the aerodynamic 
noise data available up to 1966. 

Houbolt' s curve for overall fluctuating pressure level, expressed as the 
ratio (p/q) of the mean-square fluctuating pressure level to freestream dynamic 
pressure as a function of Mach number, is shown in Fig. 1 (lower curve). 

Houbolt' s curve provides a good fit for the available data on very smooth sur- 
faces. However, investigators have measured (p/q)'s of up to 0.02 for rough- 
surfaced space launch vehicles at subsonic Mach numbers. The proposed pre- 
diction curve for p/q for attached turbulent boundary layers is shown in Fig. 1 
(upper curve). This curve was obtained by using an empirical factor of 0.02 
instead of 0.007 as used by Houbolt for the subsonic turbulent boundary layer. 

All the reliable attached turbulent boundary layer data uncovered in the literature 
search fall on or beneath this prediction curve. 

Figure 2 gives the power spectral density $(a>) in terms of p/q from Fig. 1, 
freestream velocity (U^) and boundary layer displacement thickness (6 ). This 
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curve was taken from Blake's curve (Ref. 2) for the attached turbulent boundary 
layer over a rough wall. The basic shape of the spectrum is in general agree- 
ment with the work of other investigators. 

Figures 3 and 4 give the convection velocities of the flow in the longitudinal 
direction. The broad -band convection velocities were taken from the work of 
Bull (Ref. 3) and the narrow-band convection velocities were taken from the 
empirical data of Bies (Ref. 4). 

Figures 5 and 6 give the longitudinal and lateral broad -band spatial cor- 
relations. These results of Bull (Ref. 3) are in close agreement with the work 
of a number of other investigators. 

Figures 7 and 8 give the longitudinal and lateral coherence functions. 

The coherence function is the magnitude of the cross -power spectral density 
function. Figure 7 also gives the phase angle of the cross-power spectral 
density function for the longitudinal direction. These data from Blake (Ref. 2) 
are in agreement with the work of other investigators. 

2.2 SEPARATED TURBULENT BOUNDARY LAYER 

Several investigators have studied the separated turbulent boundary layer 
associated with contour discontinuities and proturberances. Chyu and Hanley 
investigated the separated flow on an ogive cylinder in front of a 45-degree cone 
frustum. Other investigators — Speaker and Ailman (Ref. 6) and Kistler (Ref. 7) — 
studied the separated flow due to forward -facing and aft-facing steps along wind 
tunnel walls. 

Figures 9 and 10 are taken from Chyu and Hanly' s investigation. These 
figures show the mechanism of the separated flow and the general level of the 
fluctuating pressures. It is noted that the levels are higher than the levels for 
the attached turbulent boundary layer. The power spectrum of the fluctuating 
pressure of the separated turbulent boundary layer has a higher percentage of 
its energy concentrated in the lower frequencies than the power spectrum of 
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the fluctuating pressure of the attached turbulent boundary layer. The oscil- 
lating shock wave, indicated by the sharp peak in the p/q curves in Fig. 10, is 
a very localized phenomenon, with most of the power concentrated in the very 
low frequencies. 

The prediction techniques developed for the separated turbulent boundary 
layer are taken primarily from Chyu and Hanly (Ref. 5). The data associated 
with the steps were used only for enveloping the overall fluctuating pressure 
level, since these conditions are not expected on space shuttle configurations. 

Figure 11 gives the overall fluctuating pressure level as a function of 
Mach number. This curve envelopes all the data found in the literature search 
on separated turbulent boundary due to contour discontinuities and protuberances. 

Figures 12, 13 and 14 show the properties of the separated turbulent 
boundary layer based on the results of Chyu and Hanly. The power spectral 
density is shown in Fig. 12. Figure 13 presents the longitudinal coherence 
function. The lateral space correlation is shown in Fig. 14. The boundary 
layer displacement thickness indicated on these plots is the displacement thick- 
ness of the attached turbulent boundary layer immediately ahead of the separated 
region. 

2.3 BASE PRESSURE FLUCTUATIONS 

There are much less data available on base pressure fluctuations than 
there are on turbulent boundary layer. Most studies give only the overall 
fluctuating pressure level or the overall level and the power spectrum. A 
number of studies exist on the base pressure fluctuations of cone-shaped 
vehicles at low angles of attack. However, there are no data available at 
this time of the fluctuating pressure on the leeward side of a space shuttle- 
type vehicle at high angles of attack (50 to 60 degrees). Houbolt provided a 
means for determining a limiting value of the base fluctuating pressures at 
high altitudes. Houbolt theorized that since the pressure could never be 
negative, the root mean square fluctuating pressure at a point must be less 
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than 0.408 times the mean static pressure at the point. Since the static pres- 
sure at the base is generally less than the freestream static pressure, the 
maximum base fluctuating pressure can be determined as a function of alti- 
tude. Figure 15 shows this relationship. Since the space shuttle will pass 
through the high hypersonic Mach number range at altitudes over 150,000 ft, 

Fig. 15 indicates that base fluctuating pressures will not be a problem in this 
Mach number range. 

The prediction technique developed for base pressure fluctuations at 
supersonic velocities are based on the results of Robinson, Gambucci and 
George (Ref. 8), Eldred (Ref. 9), Wiley and Seidl (Ref. 10), and Lockheed 
Missiles & Space Company (Ref. 11). The data indicated that the overall 
fluctuating pressure level is a function of freestream static pressure (p^), 
instead of freestream dynamic pressure as in the case of turbulent boundary 
layer. It is proposed that a constant p/p^ of 0.05 be used for all base pres- 
sure fluctuations. This value would envelope all the data found in the literature 
search. 

Figure 16 gives the proposed prediction curve for the power spectral 
density of the base pressure fluctuations. There is an unknown constant, k^, 
associated with this curve that will have to be determined for the shuttle con- 
figuration. This results from the lack of a suitable non-dimensionalizing 
factor for the top of a shuttle-type vehicle. The spectrums presented in the 
literature were taken from bodies of revolution and used the diameter for 
the non-dimensionalizing factor. The basic shape of the spectrum, including 
rolloff rate, agrees with the results of Robinson, Gambucci and George (Ref. 8) 
and Eldred (Ref. 9). 
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Section 3 

STRUCTURAL RESPONSE 

Three general methods were initially considered for predicting the response 
of exterior panels on space shuttle-type vehicles to reentry aerodynamic excita- 
tion. They were: (1) classical dynamic equations of motion (normal mode 
approach), (2) statistical energy analysis, and (3) empirical and semi-empirical 
methods. However, two factors resulted in the elimination of empirical and 
semi-empirical methods from consideration. These factors are the complete 
lack of panel response data under reentry conditions and the complete change 
in panel structure from the Saturn vehicles to the shuttle-type vehicles that 
precluded use of Saturn test and ascent data. Consequently, solutions of the 
classical dynamic equations of motion are utilized to calculate the response in 
the low frequency modes and the statistical energy analysis is used to determine 
the high frequency response. 

3.1 LOW FREQUENCY RESPONSE 

The normal mode approach, utilizing the classical dynamic equations of 
motion, is used for the low frequency response. The method requires detailed 
information on natural frequencies, mode shapes and modal damping factors 
of the surface panels and the power spectral density and spatial correlation 
of the excitation field. 

The natural frequencies and mode shapes are calculated with the Structural 
Network Analysis Computer Program (SNAP). SNAP is a finite -element program 
developed by Lockheed -Huntsville. The structure to be analyzed is modeled as 
a network of connected beam, plate, membrane and shell elements. The program 
computes the undamped vibrational frequencies and mode shapes using an 
iterative procedure analogous to the Stodola method of beam analysis. The 
instructions for utilizing SNAP are given in Ref. 12. 
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In using the normal mode method, it is assumed that the modes are 
uncoupled, so that the total response of the structure can be determined by 
adding up the responses of the individual modes. Although this assumption is 
theoretically true only for undamped structures, it is a good approximation 
for lightly-damped structures. A problem that arises in using the normal 
mode approach is the determination of the modal damping values. They must 
be either determined experimentally or estimated. 

To derive the expression for the power spectral density of the response 
of a panel to an acoustic excitation field, we define a coordinate system for the 
panel. The following sketch shows the coordinate system. 


y 



Let x = (u, v) be the vector coordinate of a point on the panel. Let n be the 
number of modes to be utilized in the analysis. 

First, we consider the response of the panel to a single forcing function, 
F(t), applied at the point x^. Equations of motion can be set up by the use of 
Lagrange's equations. Expressed in a set of normal coordinates, 6 £ 2 * * * * * 

£ , such that the response of the system at a time t is given by 
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y(x, t) 


n 


z 

T .-1 


i r (t) * r (x) 


(i) 


Lagrange's equation gives 


£ r (t) + 2£ r 0) r 


r (t) + g r (t) 


i 

m 


L r (t) 


( 2 ) 


where 




i// r (x) = mode shape of the r mode 

= damping value for the r n mode 

th 


00^ = natural frequency of the r mode in radians per second 

2 th 

m Q (x) «/* r (x) dx = generalized mass of the r mode 


th 


L r (t) = i|< r (Xj) F(t) = generalized force of the r mode 

m Q (x) = distributed mass of the panel 

/ indicates integration over the total area of the panel. 
A 


Solving Eq. (2), we obtain 


i (t) = 

* r 


^ r (*j) F(t) 


2 2 

m (go - co + 2i £ go go) 
r r r r 


(3) 


Substituting Eq. (3) into Eq. (1), we obtain for the response at the point x^ 


y(x £ , t) = H £1 (go) F(t) 


( 4 ) 
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where 


H «<“> * ± 

1=1 


W '"r <*<> 


. 2 
m (to 
r r 


to + 2i f oo oo) 

r r 


( 5 ) 


H^(oo) is known as the transfer function. If F(t) is a random forcing function, 
the power spectral density, G (x , oo), of the displacement, y(x ), is given by 

y * t 


^y( X ^, ^0) 


H fl (to) 


G f (oo) 


( 6 ) 


where 


G f (oo) = power spectral density of the forcing function F(t). 

Next, we consider a number of forces, Fj(t), F 2 (t), . .., F^(t), acting 
at the points, Xj, x-^, . . . x^, respectively. In this case, the relation shown 
in Eq. (6) becomes more complex due to possible correlations between the 
forcing functions as well as the increased number of forcing functions. For 
this case, the power spectral density of the response at x^ becomes 

= £ £ H *ta «“> <“> G ojJ «“> (7 > 

0=1 P =1 


where 


G 


«P 


(M) 


cross -power spectral density of 


F and 
a 



* 


Indicates the complex conjugate. 
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We now return to the case of interest, where the forcing function is an 
acoustic excitation field. For this case, the double summation of Eq. (7) be- 
comes a double integral to give 


G y (x C M) = / / H ia < M > H £3 <“> G ap (M) dx a dx 3 • 
A A 


th 


Modal impedance, Z r (CD), of the r mode is defined by 


2 2 

Z (CD) = m (CD - CD + 2i t CD CD) . 

r r r *r r 


(8) 


(9) 


Combining Eqs. (5), (8) and (9) and simplifying gives 


G y( X jf, 



n ^(x,) 

s=l Z r ( '- 0) Z s <“> 



(cd) dx 

a 



(10) 


Separating Eq. (10) into equal-mode terms and unequal-mode terms and multi - 
2 

plying by ^ ^ , where 4 >(cd) is the power spectral density of the acoustic 

A d 4>(cd) 

excitation field, gives 


n 


G (x^CD) = A 4 >(cd) ^ 

r= 1 


*1 < x <> 


Z r (&>) 


// 

A A 


iff (x ) \l/ (x ) G 
r r v a r r v 3 7 

.2 


$(cd) a 3 


(H) 


2 A A ^r (x P f /* ’M X a ) ^s (x R ) G oR (a)) - - 

+ A^ 4>(cd) V V T / / — A s . J dx dx . 

£ .4i z>) z,(“) 4 a 6 

s^r 
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The joint acceptance squared, j (co), is a measure of the correlation between 
the pressure field and the mode deflection. It is given by 


.2 

J r 


(CO) 


ff W»r<*B> 

• A J A AT 


G qB |M> 
4>(co) 




(12) 


The cross -joint acceptance squared, j (co),is defined by 

IT S 




( 60 ) 




(13) 


Substituting Eqs. (12) and (13) into Eq. (11) gives 


Gy(x^, co) = A 4>(co) 


^( X f) Jr {a>) 

2 


r=l 


Z r (C0) 


2 JL JL ^ r ( x P (*#) j^s (£i)) 

+ A Z 4»(co) V V r £ S £ rS 

l l z (60) Z (CO) 

r=l s=l r s' 

s^r 


(14) 


Since the cross-joint acceptance is generally very small compared with the 
joint acceptance, the second term in Eq. (14) is dropped to give 


G y ( x f » ^) 


A 6 4>(co) ^ 

r=l 


^r jr 




(15) 


Although the cross-power spectral density function, (CO), in Eq. (12) is a 

complex function, the imaginary part integrates to zero. Therefore, G^ (co) 
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can be replaced by its real component. But the real component of G (03) 
divided by ^(co) defines the spatial correlation function, R^ (go). Therefore, 
Eq. (12) simplifies to 


4 <®> 



ilj (x ) (x. ) R . (03) dx dx . 

r r a ^r (3 a(3 a 0 


(16) 


Taking the absolute value of Eq. (9) and squaring gives 


Z (03) 
r 


2 2 
= m 

r 


[«■£ " ^ 2 ) 2 + 03 2 ] . 


(17) 


The power spectral density of the acceleration response, G.. (x^, 0)),is given by 


Gy (x^, - G q 03 Gy ( X £» 


(18) 


where C is a conversion constant that depends on the units of the power spectral 
o 2 2 

densities. If G is in inches /radian and Gy is desired in g /radian, then 

C = 0.6709 x 10“ 5 . 
o 


A Fortran computer program has been written to calculate the power 
spectral density of the displacement response and the acceleration response 
of a panel, utilizing the output tape from SNAP for the natural frequencies and 
mode shapes. The program is based on Eqs. (15) through (18), with the double 
integration being approximated by a double summation. A listing of the program 
is given in Appendix B. 
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3.2 HIGH FREQUENCY RESPONSE 


Statistical energy analysis is used to estimate the high frequency response 
of the exterior panels on space shuttle-type vehicles. The statistical energy 
approach to response estimation was developed to handle high frequency vibra- 
tions. The method basically works with averages and eliminates the need for 
detailed information on natural frequencies and mode shapes. Lyon and 
Maidanik first postulated the technique in Ref. 13. Further work in this field 
was done by a number of other investigators, including Smith and Lyon (Ref. 

14), Scharton (Ref. 15), and Ungar (Ref. 16). The following methods developed 
for space-shuttle panels are based on Ungar' s results. 


We consider a panel exposed to a diffuse sound field in a reverberation 

chamber. The total kinetic energy, K (to., 6.), in the chamber in a frequency 

a J J 

band, 6., centered at frequency to. is given by 
J J 


p (CO., 6 .) V 

K (0> ,5.) = L 

J J 2 n c* 


(19) 


where 


c = speed of sound in the acoustic medium 


p (to., 6.) = mean- square fluctuating pressure in the frequency 
** band, 6^, centered at frequency co^ 

V = volume of the acoustic medium 


p a = density of the acoustic medium. 


The average kinetic energy per acoustic mode, <^K (to., S.)> , in the frequency 

a j j m 


band, 6., centered at frequency 


60^ is given by 
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< K (co., 6 .)> 
a' y y'm. 


K (co„ 6.) 

a j i 

n (co.) 6. 
a' r J 


( 20 ) 


where 


n (co.) = modal density of the acoustic medium at frequency co.. 
a J J 


The modal density, n (co.), of the acoustic medium is given by 

a J 


n 

a 



V u? 

.-rV 

2?r c 


( 21 ) 


Substituting Eqs. (19) and (21) into Eq. (20) and simplifying gives 


2 2 . . . 

it c p (co , 6.) 

<K (co., 6 .)> = = — J — J- 

a' j’ j"m 


2 „ 
p co. 6. 
M a j j 


( 22 ) 


The total kinetic energy, K (co., 6.), of the panel in a frequency band, 

P J J 

6^, centered at frequency ox is given by 


m a^(co., 6.) 
K (co., 6.) = — E — 5 - J — J- 

P J J 

J 


(23) 


where 


a 2 (co., 6.) 
J J 

m 

P 


mean- square acceleration of the panel in the frequency 
band, 6^, centered at frequency C0j. 

mass of the panel. 
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The average kinetic energy per panel mode, <K (co., 6.)> , in the frequency 

P 3 3 m 

band, 6., centered at frequency co. is given by 
3 3 


K (CO., 6.) 

<K (co., 6.)> = — P . X - J - 

P 3 3 m n p j 6 j 


(24) 


where 


n (co.) = modal density of the panel at frequency CO.. 
P 3 3 

Substituting Eq. (23) into Eq. (24) gives 


<K (CO., 6 .)> 
p j’ j'-'m 


_2 

m a (co., 6.) 

e I_.^L 

2 co. n (CO.) 6. 

3 P 3 3 


(25) 


Ungar (Ref. 16) has shown that for two coupled sets of modes where energy 
is supplied to one mode set and no energy is supplied to the other set, the ratio 
of the average modal kinetic energy of the two sets can be expressed in terms of 
a coupling loss factor and a modal loss factor. Applying Ungar 1 s formulation 
to the present problem gives 


<K (co . , 6 . )>■ 

P v i l m 

<K (C0 , 6 )> 
x a j’ j 'm 


ri (co., 6.) 
. P . a . J J L 


■n (co., 6.) + ri (CO., 5.) 
P a 3 3 ! P 3 3 


(26) 


where 


r) (co., 6.) = average coupling loss factor for energy flow from 
P a ** a panel mode to an acoustic mode in the frequency 
interval, 6 ., centered at frequency CO. 
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*1 


(0)., 6.) 
P J ) 


average modal loss factor for a panel mode in the 

frequency interval, 6., centered at frequency 0).. 

J J 


Substituting Eqs. (22) and (25) into Eq. (26) and simplifying gives 


-2 

a 


( 0 )., 6 .) 

' J J 



2 

c n (CO. ) q (CO., 6.) p (CO., 6.) 

P .1 Pa .1 J J -J 

p [q (OX, 6.) + q (CO., 6.) 

a L 'pa' j J P J J J 


(27) 


Equation (27) provides a method for determining the acceleration response 
of a panel excited by a reverberant acoustic field. To apply this formula, one 
must be able to divide the frequency range of interest for the response into a 
number of frequency intervals. These intervals must be chosen so that the 
structural modes in a given interval have similar modal loss factors. All 
structural modes in the given interval must have equivalent coupling to all 
acoustic modes in the interval, and the coupling loss factors must be similar. 
For each frequency interval, one must have a good estimate of the modal 
density of the panel, the average modal loss factor of the panel modes, and the 
average coupling loss factor for energy flow from a panel mode to an acoustic 
mode. 
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Section 4 

RESPONSE OF TPS TEST PANEL 

A test panel of a proposed thermal protective system (TPS) configuration 
was fabricated at Marshall Space Flight Center for reverberation chamber 
tests. The test panel, which is made of Haynes 25 Alloy L-605, is shown in 
Fig. 17. It consists of two individual corrugated panels which overlap each 
other by one inch. Only one half of one of the panels is shown in Fig. 17. 

The panel is symmetrical about the A -A axis. The panel is attached to a 
corrugated stand off -panel along one edge and a series of flexible clips along 
the other edge. Only half of each clip is shown. The clips are symmetrical 
about the B-B axis, with the other panel attached to the other half of each 
clip. Each panel has a series of horizontal stiffeners parallel to the corruga- 
tions. The stiffeners are attached to the panel by a series of spot welds. The 
panel and stiffeners are attached to the standoff -panel and clips by bolts (not 
shown). The standoff -panels and clips are bolted to large aluminum channels 
perpendicular to the corrugations. 

To determine the low frequency response of this panel in the reverbera- 
tion chamber, a finite element model of the panel was developed. In the model, 
the panel is represented as a series of plate elements and the stand off -panel, 
clips and stiffeners are represented as a series of beam elements. The points 
of connection between the aluminum channels and the standoff -panel and clips 
are considered fixed in the model. Since the two individual panels are con- 
nected at fixed points on the clips and by a sliding overlap, they can be con- 
sidered decoupled. Therefore, it is sufficient to model only one of the panels. 
A plane -of- symmetry technique is used so that only one-half of each individual 
panel is modeled. The original model is shown in Figs. 18, 19 and 20. 

The model was analyzed by the Structural Network Analysis Program 
(SNAP) to determine the natural frequencies and mode shapes of the low 
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frequency modes. However, four out of the first five modes turned out to be 
local modes of one cantilever overhang of the panel. The first mode is shown 
in Figs. 21, 22 and 23. Since the local cantilever modes do not appreciably 
affect the response near the center of the panel, the model was changed 
to eliminate the overhangs. The revised model is shown in Figs. 24, 25 
and 26. 

The revised model was analyzed by SNAP to determine the first five 
modes. All five modes are panel modes. They are shown in Figs. 27 through 
41. 


A computer program was written to compute the response at any designated 
point on the panel. The program uses the output tape from SNAP for the natural 
frequencies and mode shapes. The spatial correlation used for the reverberant 
field is given by 



(CO) 


sin(A 





co/c) 

co/c 


(28) 


where 



distance between points a and (3 . 


A description of the program is given in Appendix B. 


The acoustic test criteria for the reverberation chamber in one-third 
octave bands is given in the following table. The power spectral density (PSD) 
of the acoustic test criteria was determined from the one-third octave band 
data. The PSD is shown in Fig. 42 on a log -log scale. 
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TABULATION OF ACOUSTIC TEST CRITERIA 


One -Third Octave Band Acoustical Specification in dB re 2 x 10 


n/m 2 


Geometric 

Mean Freq. 

dB 

(Hz) 

5.0 

135.0 

6.3 

136.0 

8.0 

137.0 

10.0 

137.5 

12.5 

138.5 

16.0 

139.5 

20.0 

140.0 

25.0 

140.5 

31.5 

141.0 

40.0 

141.5 

50.0 

142.5 

63.0 

144.5 

80.0 

146.5 

100.0 

147.5 

125.0 

148.5 

160.0 

148.5 

200.0 

147.5 

250.0 

146.5 

315.0 

145.5 

400.0 

144.0 

500.0 

143.0 

630.0 

142.5 

800.0 

142.5 

1000.0 

141.5 

1250.0 

141.0 

1600.0 

140.5 

2000.0 

139.5 

2500.0 

138.5 

3150.0 

137.5 

4000.0 

137.0 

5000.0 

136.0 

6300.0 

135.0 

8000.0 

134.0 

10000.0 

133.0 


Overall SPL 


158.0 
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The response program was utilized to calculate the response of the panel 
at the point P (Fig. 17), which is over a stiffener near the center of the panel. 

The revised model of the panel (without overhangs) was used. Two separate 
runs were made, one using 4% damping for all modes and one using 5% damping 
for all modes. The results obtained with 4% damping are shown in Figs. 43 
through 46. The results of the 5% damping run are given in Figs. 47 through 
50. These figures give the power spectral densities of the response displace- 
ment and acceleration using both linear and log -log scales. The overall variance 
and the root mean square (RMS) of the response is also given on each figure. 

The response power spectral density is dominated by the response of the 
second mode, at approximately 198 Hz. Analysis of the mode shapes of the 
panel and the spatial correlation of the acoustic field reveals that the joint 
acceptance squared is much higher for the second mode than for any of the other 
modes. This effect overshadows other factors, such as the higher acoustic 
excitation level at the natural frequency of the first mode. 

Comparison of the results of the two runs shows the effect of damping. 

Due to the predominance of the second mode, the damping value of the second 
mode is the primary determining factor for the height of the resonant peak 
and the RMS of the response. 
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Fig. 1 - Overall Fluctuating Pressure Level for 
Attached Turbulent Boundary Layer 
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Fig. 6 - Lateral Broad-Band Spatial Correlations for 

Attached Turbulent Boundary Layer (from Bull, 

Ref. 3) 
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Fig. 7 - Longitudinal Coherence Function for Attached Turbulent 
Boundary Layer (from Blake, Ref. 2) 
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Fig. 9 - Separated Turbulent Flow Ahead of Cone Frustum 
(from Chyu & Hanly, Ref. 5) 



Fig. 10- Longitudinal Variation of Fluctuating Pressure in Separated Turbulent 
Flow (from Chyu & Hanly, Ref. 5) 
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Fig. 15 -Maximum Base Fluctuating Pressure by Altitude 
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Fig. 18 - Original Model of TPS Test Panel, View 1 
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Fig. 19 - Original Model of TPS Test Panel, View 2 
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Fig. 21 - First Mode of the Original Model, View 1 
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Fig. 22 - First Mode of the Original Model, View 2 
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Fig. 23 - First Mode of the Original Model, View 3 
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Fig. 24 - Revised Model of TPS Test Panel, View 1 


54 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225410 


UNDEFORMED STRUCTURE 



pi.O r i 
2 


T p S TEST PANEL * L-605 HAVNES 25 


e 


SC A 


Fig. 25 - Revised Model of TPS Test Panel, View 2 
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Fig. 27 - First Mode of the Revised Model, View 1 
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Fig. 28 - First Mode of the Revised Model, View 2 
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Fig. 29 - First Mode of the Revised Model, View 3 
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Fig. 30 - Second Mode of the Revised Model, View 1 
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Fig. 31 - Second Mode of the Revised Model, View 2 
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Fig. 32 - Second Mode of the Revised Model, View 3 
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Fig. 33 - Third Mode of the Revised Model, View 1 
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Fig. 34 - Third Mode of the Revised Model, View 2 
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Fig. 36 - Fourth Mode of the Revised Model, View 1 
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Fig. 37 - Fourth Mode of the Revised Model, View 2 
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Fig. 38 - Fourth Mode of the Revised Model, View 3 
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Fig. 39 - Fifth Mode of the Revised Model, View 1 
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Fig. 40 - Fifth Mode of the Revised Model, View 2 
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TPS PANEL HAYNES 25 ALLOY L-605 VAR = 0.74332 E-04 

DISPLACEMENT RESPONSE IN INCHES JOINT 53 RMS « 0.86216 E-02 



Fig. 43 - Power Spectral Density of the Response Displacement at P, 

Linear Scale, 4% Modal Damping 
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TPS PANEL HAYNES 25 ALLOY L-605 VAR * 0.74332 E-04 

DISPLACEMENT RESPONSE IN INCHES JOINT 53 RMS * 0.86216 E-02 



Fig. 44 - Power Spectral Density of the Response Displacement at P, 
Log -Log Scale, 4% Modal Damping 
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ACCELERATION RESPONSE IN S*S JOINT 53 R«S s 0.32164 E 02 



Fig. 45 - Power Spectral Density of the Response Acceleration at P, 

Linear Scale, 4% Modal Damping 
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ACCELERATION RESPONSE IN 6*S JOINT 53 RHS = 0.32164 E 02 



Fig. 46 - Power Spectral Density of the Response Acceleration at P, 
Log -Log Scale, 4% Modal Damping 
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TPS PANEL HAYNES 25 ALLOY L-605 VAR = 0.60338 E-04 

DISPLACEMENT RESPONSE IN INCHES JOINT 53 RMS = 0.77678 E-02 



Fig. 47 - Power Spectral Density of the Response Displacement at P, 
Linear Scale, 5% Modal Damping 
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TPS PANEL HAYNES 25 ALLOY L-605 VAR = 0.60338 E-04 

DISPLACEMENT RESPONSE IN INCHES JOINT 53 RMS = 0.77678 E-02 



Fig. 48 - Power Spectral Density of the Response Displacement at P, 
Log -Log Scale, 5% Modal Damping 
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Fig. 49 - Power Spectral Density of the Response Acceleration at P, 

Linear Scale, 5% Modal Damping 
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Fig. 50 - Power Spectral Density of the Response Acceleration at P, 

Log -Log Scale, 5% Modal Damping 
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Appendix A 


A. 1 INTRODUCTION 

A large portion of the surface of space shuttle -type vehicles will be 
subjected to attached turbulent boundary layer during reentry into the earth's 
atmosphere. An investigation was made of the present state of knowledge 
of the turbulent boundary layer. A theoretical description of the boundary 
layer is presented. Problems encountered in measuring turbulent boundary 
layer parameters are discussed, along with results of experimental investiga- 
tions. 

A. 2 DISCUSSION 

For a perfect fluid, i.e., frictionless and incompressible, motion does 
not create any tangential forces between contacting layers of fluid. Because 
of the absence of these forces, a slip condition must be assumed between the 
moving fluid and a solid boundary, i.e., a difference in relative velocities. 

For real fluids, however, tangential or shear forces can be transmitted be- 
tween layers due to the viscosity. Then, because of the existence of intermo- 
lecular attractions, the fluid will adhere to a solid wall and give rise to a 
shearing stress which is transmitted through the fluid. In a thin region 
near the solid boundary the velocity increases from zero to its full value 
which corresponds to the external frictionless flow. The thickness of this 
layer increases along the boundary in the downstream direction and increases 
with an increase in viscosity (v). Initially the particles of fluid move with a 
uniform velocity along straight paths. The flow is well ordered and the 
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moving particles can be visualized to form laminae. As the Reynolds number 
(V^jf/t/)is increased the orderly pattern ceases and strong mixing of the par- 
ticles occurs. The pattern of streamlines at a fixed point, as well as the 
velocity and pre ssure, becomes subjected to continuous fluctuations and there 
exists a continuous transport of energy from the main flow into the large eddies. 
Energy is dissipated, however, by the small eddies. This dissipation process 
occurs in a narrow strip inside the boundary layer in the region near the solid 
boundary. The eddies or balls of fluid that are formed have their own intrinsic 
motion which is superimposed on the main flow. The size of the eddies which 
are continually formed and dissipated determines the scale of the turbulence, 
i.e., the size is determined by the external condition of the flow. 

° As would be expected, the velocity at which the eddies are conveeted 
downstream varies through the boundary layer. Evidence of this is presented 
in Schlichting (Ref. A.l), from Nikuradse (Ref. A. 2) and Tollmien (Ref. A. 3), in 
photographs taken with a camera which moved with the flow at varying speeds. 
From correlation of velocity and pressure measurements Willmarth and 
Wooldridge (Ref. A. 4), Bull (Ref. A. 5), Serafini (Ref. A. 6) and Blake (Ref. A. 7) have 
shown that the convective velocity of the eddies varies with position in the bound- 
ary layer and can be associated with the eddy size. These data were obtained 
from measurements in the boundary layer on a flat plate with zero pressure 
gradient. The average convective velocity obtained from the broadband space- 
time correlations of a number of pressure transducers is presented in Fig. A-l 
as a function of the separation distance (r) between transducers. For small 
values of r the convection velocity ratio U c /U begins to decrease significantly 
as r decreases. These results can be reasonably interpreted by considering 
that the smaller eddies have a shorter lifetime. Then as the distance between 
transducers increases, the correlation of the pressure fluctuations decreases, 
since these are only affected by the larger scale eddies. This implies that the 
small eddies are concentrated closer to the wall, since these are conveeted at 
slower velocities. The eddies in the boundary layer can then be visualized in 
two limiting regions; one near the wall which contains the smaller scale eddies 
that travel at low speeds and the other near the outer region which contain the 
larger scale eddies and move with a velocity approaching U . 
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The physical description of the boundary layer is helpful in considering the 
results and problems that would occur in measuring the unsteady pressure. 
The sketch below shows a pressure transducer flush mounted to a flat plate 
immersed in a viscous flow. Typical nomenclature definitions are presented 

below. 



d = diameter of the eddy 
e 

D = diameter of the transducer 
u = circular frequency 

jJ{ 

6 = displacement thickness of the boundary layer 

♦ 

k = u6 /U (reduced frequency) 

U = convective eddy velocity 
c 

The pressure distribution in the eddy should be similar to the Rankine 
combined vortex shown on the following page. 
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r /R 
e' 


Therefore, the velocity and diameter of the eddy as it moves across the pres- 
sure transducer should be related to the frequency of the output as 

U 

c 

W = d“ 
e 

If the transducer diameter is larger than one-half the wave length, the meas- 
urement will be attenuated. 

The energy at various frequencies can also be considered qualitatively 

* , 

by examining the expression for the reduced frequency(k = u6 /U ), 

* 

assuming 6 = 0(10 6 ) 


and letting 


then 


(J 



U 

oo 


e 




k 


if k = 0.1 then K x =0(d e /6). However, as k decreases, i.e., if k « 1 then 
d /6 must be larger than 1 for reasonable values of K^. This implies Kj 
approaches 1. Therefore, setting K 1 = 1 
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d 

e 

5 


1_ 

= 10k 


and for k « 1 implies d g > 6. 

The probability of eddies with diameters larger than 6 is relatively 
small and the energy measured by the transducer at these reduced fre- 
quencies should be much lower. 

Also, k » 1 implies that d g is very small for reasonable values of 
(i.e., Kj < 1), and the energy content of these eddies sensed by the transducer 
should be significantly lower. Therefore, from this model, the PSD of the 
measured pressure should be a maximum between 0.1 < k < 1.0 and orders 
of magnitude lower when k > 10 and k < 0.01. 

The above examination of the physical description also aids in defining 
the similarity parameters for collapsing the available data. The above analy- 
sis suggests that the boundary layer thickness which is related to the eddy 
diameter (the displacement (6 ) or momentum (6g) thickness) should be used 

to nondimensionalize length. This leads to a reduced frequency defined as 

♦ 

k = u6 /U where U is the local velocity of the external frictionless flow. Ex- 
perimental and theoretical information on the pressure levels shows reason- 
ably good correlation has been found using P rms / ( 1 as a parameter. Lilley 
(Ref. A. 8) suggests that P rms / T w wou ld be a suitable theoretical nondimension- 
alized parameter, where t is the wall shear stress. This theoretical pre- 
diction is dependent on the similarity of the boundary layer velocity profiles. 
Neither the effects on pressure fluctuations of roughness (Refs. A. 7 and A. 9) nor 
those of pressure gradient (Ref. A. 10) can be accurately predicted by using t , 
although some improvement can be seen. Thus, use of T increases the 
complexity of prediction without any significant advantages. 

A number of both wind tunnel and flight test programs have been con- 
ducted to measure the fluctuating pressures beneath a turbulent boundary 
layer for subsonic flow. These data are inconsistent, however, and a number 
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of difficulties are presented in evaluating the cause of the anomalies. There 
is little or no uniformity in nondimensionalizing the data and in presenting 
the results. In most cases of the reported flight measurements no informa- 
tion is furnished about the boundary layer or local flow parameters. These 
data are normally nondimen sionalized by calculated values of the boundary 
layer thickness and the freestream parameters independent of the body loca- 
tion. Flight measurements are also subject to background noise caused by the 
engines and turbulence generated by sections of the aircraft upstream of the 
measurement location. Although the flow conditions are well described for 
the wind tunnel measurements, the tunnel noise and turbulence can have large 
effects on the data and both flight and wind tunnel results are subject to errors 
caused by the frequency limitations of the instrumentation and the inherent 
loss in sensitivity of a finite size pressure transducer. A summary of the 

measured values of P /q are presented in Table A-l. 

rms 4 5 * 7 oo 

After reviewing these reports and observing the scatter in the data, 
calculating an average value or performing statistical treatment of all the re- 
ported values would be meaningless. The scatter in these data far exceed 
the possible errors reported by any of the authors. The scatter in the data 
is not significantly reduced, even upon a selected compilation where data 
suspected to have attenuation of the high frequency components caused by the 
data acquisition equipment or transducer size are discarded. It was therefore 
concluded that the large scatter was either due to error or caused by a subtle 
difference in the testing techniques or conditions. 

The results of the various investigations were further studied only if 
the following conditions were met 

1. D/6* < 0.5, 

2. Data acquisition equipment with the capacity to accurately 

record results for large k = U)6*/U was used, 

00 

3. Low background noise, 

4. Known (measured) values of the local flow conditions were 
available, and 

5. Details of experimental procedures were reported. 
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Table A- 


SUMMARY OF EXPERI 


R e po r t 

P /q 

rms' oc 

M 

oo 

3P/3S 

(psi/in.) 

V 

(Range; 

o 

N 

°V 

Speaker 
*> Ailman 
(Ref. A. 1,1) 

0.0058 

0.0059 

0.0070 

0.0075 

0.0044 

0.42 

-9x 10‘ 4 

m 

3.40 

0.66 

0.66 

0.60 

8.00 

0.0052 

0.0057 

0.0050 

0.0070 

0.0067 

0.59 

-9x 10' 4 



Serafini 
(Ref. A. 6) 

0.0074 

0.0072 

0.0077 

0.0068 

0.0058 

0.0054 

0.0048 

0.0045 

0.58 

-1.1 x 10" 3 

El||^ 

1 

0.20 

0.20 

0.20 

0.30 

1.2 

3.4 

5.5 
9.0 

Bull 

(Ref. A . 5) 

0.0045 

to 

0.0053 

0.3 

0 

1 to 
6x10^ 

0.3 

to 

1.0 

W ilby, Bhat 
& Gloyna 
(Ref. A .:"'l 2) 

0.0053 

0.78 


lxlO 5 

a 

Schloomer 
(Rof.A. 10) 

0.0052 

0.0046 

0.07 to 

0.09 

0 

10 4 

0.406 

0.0075 

0-00585 

0.09 to 
0.13 

"4.4 x 10' 3 

1.8 x 10 4 

0.292 

0.00362 

0.00502 

0.14 

0.14 

BH 

Hi 

2.36 

Ludwig 
(Ref. A'. 13) 

0.008 

to 

0.011 

0.17 

to 

0.02 



0.092 

to 

0.74 

Gibson 
{Ref. A . 14) 

0.0035 

to 

0.01 

0.24 

to 

0.80 



0.36 

to 

1.54 

Leech 
(Ref. A . 15) 

0.003 

to 

0.01 

0.90 

to 

1.16 


B9 

■ 

Willmarth & 
W oold ridge 
(Ref. A . 4) 

0.0056 

0.0052 

0.0049 

0.0047 

0.19 

0 

4 to 

5.2 xlO 4 

0.12 

0.21 

0.33 

0.45 

W iley & 

Seid 1 

(Ref. A . 16) 

0.006 

to 

0.02 
Sub- 
sonic 
V alues 
- 0.01 

■ 
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i 

MENTAL DATA 


Corrected 
for Trans- 
ducer Size 



'u_6 

U 


Yes 


7x 10' 7 


Comments 


= 5 


Wind Tunnel Test 


Yes 

3 to 5 x 10" 6 


No 

2x 10' 6 

Wind Tunnel Test 

No 

8 x 10' 7 

Wind Tunnel Test 

Yes 

5x 10' 7 

Flight Test 

Yes 

No 

6x 10' 7 

Wind Tunnel Test 

High F requency 

Cutoff at 10 kHz 
♦ 

or u6 /U = 8 and 
00 * 

Low Cutoff u>6 /U ^ 0.2 
00 

Yes 

No 

6 x 10" 7 

Low Frequency 

Cutoff at u>6*/U =0.12 

00 

No 

Yes 

6x 10' 7 

Transducer Very 
Large 

No 

1 to8x 10* 4 

Channel or Duct Flow 

No 


Local Condition Un- 
known (Flight Test) 
Large Scatter in Data 

No 


Flight Test (F 102A 
Aircraft) Frequency 
Limited to 2 kHz 

No 

1 x 10' 6 

Wind Tunnel Test 
Low Frequency Cutoff 

u6*/U =0.13 (High Low 

00 

Frequency Noise Level) 




Wind Tunnel Test of 
X 20 Vehicle 
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Results of those studies meeting most of the above conditions are tabu- 


lated below. 


Report 

Mach 

No. 

P /q 

rms' oo 

R 

e 6* 

(Range) 

/ * 
D/6 

Corrected 
for Trans- 
ducer Size 

Willmarth 
(Ref. 4) 

0.19 

0.0056 

0.0052 

0.0049 

0.0047 

0.0050 

4 to 5 .2 x 10 4 

m 

No 

Bull 
(Ref. 5) 


0.0045 to 
0.0053 

1 to 6 x 1 0^ 

1 .0 to 
0.3 

No 

Blake 
(Ref. 7) 

0.07 

0.15 

0.00876 

1 to 2 x 10^ 

0.11 

No 

Speaker 
(Ref. 11) 

0.42 

0.59 

0.0044 to 
0.0075 

2.1 to . 
2.4 x 10 

8.0 to 
0.60 

Yes 

Serafini 
(Ref. 6) 

0.58 

0.0045 to 
0.0075 

2.5 to 
8 x 10 

9.0 to 
0.2 

No 


NOTE: 9P/3S % 0 for all the above. 


9|C 

The wall pressure spectrum results with the exception of Blake are presented 
in Fig. A-2. 


As seen in the above table, there is much scatter in the P / q values. 

rms' oo 

Some of the scatter can be attributed to transducer size; however, if values of 

P /q are compared with results of Willmarth, Bull, Blake and Serafini 
rms' ^oo _ 

for D/6 of 0.1 to 0.2, ^ > rms / < l 00 still ranges from 0.005 to 0.00876. The only 
plausible explanation that could be found is the effect of surface roughness. 
Willmarth and Wooldridge (Ref. A. 4) flush-mounted the transducers to a one- 
inch thick (oil-lapped) steel plate 20 inches in diameter, which was supported 
independently of the tunnel by a large pedestal. The instrumented plate was 
flush with the tunnel floor which was a newly installed, varnished and waxed 


j|( 

The_pressure spectrum presented by Blake was inconsistent with his results 

of P /q . Blake also presented the spectrum obtained by Bull and Willmarth 
rms' oo 

in the same figure and it appears as though all the data were plotted incorrectly. 
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sheet of masonite extending 14 feet upstream. The holes in the plate for 

mounting the transducers, which were not in use, were filled with brass plugs. 

These were fitted within _+0. 001 inches of the surface. On a second series of 

tests in which even greater care was taken to ensure smooth flush surfaces, 

a reduction of approximately 13% in the measured value of P /q was noted. 

Bull (Ref. A. 5) reported that the tunnel walls were ground and polished and the 

transducers mounted in six-inch diameter plugs which were attached to the 

wall so that the transducer and wall formed a continuous surface. The values 

of p rrr .„AL measured by both authors were nearly the same (P /q % 0.0053 
•this oo 7 rms' oo 

jf 0.0003). Serafini (Ref. A. 6) does not describe the details of the tunnel nor the 
installation of the transducers other than that they are flush-mounted. He im- 
plies, however, that the tunnel walls were smooth, but not polished. Speaker 
and Ailman, also, do not describe the installation; however, from photographs 
presented, the transducers appear to be flush-mounted on an unpolished alumi- 
num plate which had a number of fine scratches on the surface. 

Blake (Ref. A. 7) flush-mounted the transducers in an aluminum traverse 
which was flush-mounted with the test surface consisting of a sheet of formica- 
coated plywood. The aluminum traverse disc was drilled along a diameter for 
positioning the microphones at various locations. The holes were filled when 
the disc was not in use; the assembly was reportedly smooth to the touch. A 
more exact description of the conditions of the surface was not given. Blake 
also investigated the effect of surface roughness by adhering sand to the test 

surface and the transducer diaphragm. These results are presented in the 
table below. 

Wall Designation k^ 

Smooth 

S-S 0.0563" 

D-S 0.0563" 

D-L 0.092" 


k 

s 

P /q 

rms' oo 

P /■ 

rms' 


0.00876 

3.59 

0.106" 

0.0202 

3.83 

0.0736" 

0.01605 

3.2 

0.192" 

0.0184 

2.87 
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NOTE: S-S = Small, sparsely packed sand particles 

D-S = Small, densely packed sand particles 

D-L = Large, densely packed sand particles 

k = Mean geometric roughness height 
g 

k = Equivalent hydrodynamic sand roughness height 


Coe (Ref. A. 17) shows large changes in the pressure fluctuations measured by 
a transducer as it was moved in and out from a flush position and pressures 
measured by a flush-mounted transducer positioned 0.118 6 downstream. For 
a change in level of the variable transducer of 0.002 y/6, the output changed 
approximately 60 to 70%, while the output of the fixed transducer changed 
approximately 10%. Results of Speaker and Ailman are again of interest 
since a large number of transducers were aligned in the direction of the flow 
and the surface was not polished. The P values are presented in 

Figs. A-3 and A-4 as a function of the downstream location. Although it can be 
argued that the increase in P /q with distance is due to the change in 
D/6 , the possibility of the progressive effects of the relatively rough surface 
and transducer misalignments also exist. 

The narrow-band convection velocities as measured by Blake (Ref. A. 7) 

jjC . 

and Bull (Ref. A. 5) are presented in Fig. A-5. They agree well for <J6 /U^ > 0.3 

however, for oj 6*/U < 0.3, the pressure fluctuations have low phase velocities 

00 

which are strongly dependent on the microphone separation distance used in 
their measurement. The association of these pressure disturbances with a 
particular eddy system within the boundary layer is difficult. Also, from the 
earlier analysis of the equation Kj = 10 d^k/6, for k =0(0.1), both high- and 
low- velocity eddies could contribute to the pressure fluctuations, and large 
scatter in the measured narrowband convective velocities might be expected. 
However, it was also shown that as k decreases, the contribution should come 
from the large, high-speed eddies. 


A-10 


LOCKHEED- HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225410 


A. 3 CONCLUSIONS 


Results of the experimental investigations of the characteristics of the 
fluctuating pressure measured beneath a subsonic turbulent boundary layer 
are in relatively good agreement for measurements of the broad and narrow 
band convective velocities; unfortunately, however, in measurements of the 
spectrum and rms levels of the pressure, the data were scattered. Most of 
the scatter can be attributed to the attenuation or error caused by the finite 
size of the transducer or high frequency limitations of the data acquisition 
systems. However, there is still considerable scatter for data which were 
obtained when small and approximately equal transducers (D/6 = 0.1 to 0.2) 
were used and which were recorded by equipment capable of obtaining the 
high-frequency components. The data obtained by Blake (Ref. A. 7) show large 
effects due to surface roughness and data by Coe (Ref. A. 17) demonstrate large 
effects caused by poor alignment of the transducer. This, combined with the 
variation in the smoothness of the surfaces used in the various investigations, 
indicates that the quality of the surface could have caused the observed scatter 
in the data. 
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Fig. A - 2 - Mean-Square Spectra of Wall Pressure Fluctuations 
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O Unthickened Boundary Layer) Piezoelectric 
A Thickened Boundary Layer j Transducer D = 0.06 
Filled Symbols — Condenser Microphone (FM) D = 0,25 
Tagged Symbols — Condenser Microphone (AM) D = 0.25 



Distance Along Tunnel Wall, X (in.) 
Flow w 


Fig. A -3 - Variation of the rms Pressure Fluctuations Along the Tunnel Wall, M = 0.42 
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Unthickened Boundary Layer/ Piezoelectric 

( Transducer D = 
Thickened Boundary Layer ) 

Filled Symbols — Condenser Microphone (FM) D 
Tagged Symbols — Condenser Microphone (AM) D 


Flow 


Distance Along Tunnel Wall, X (in.) 
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Fig. A -4 - Variation of the rms Pressure Fluctuations Along the Tunnel Wall, M 
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Appendix B 


B. 1 DISCUSSION 

This program is written in Fortran V for the Univac 1108 computer at 
Marshall Space Flight Center. It uses as input the output tape from the Struc- 
tural Network Analysis Program (SNAP) along with data cards. The SNAP 
output tape provides the natural frequencies and mode shapes, normalized to 
unit generalized mass, of a panel structure to be exposed to a reverberant 
acoustic field. The data cards provide the power spectral density of the acoustic 
field, the damping values of the panel modes, designation of the panel modes 
to be used in the analysis, designation of the joints in the finite element model 
to be used in the double area integration over the mode shapes, and designation 
of the joint where the response is to be determined. The primary output of 
the program is four SC 4020 plots, giving the power spectral densities of the 
response displacement and acceleration using both linear and log-log scales. 
Although the program was written for a reverberant acoustic field, it can be 
easily modified to handle other types of acoustic fields by changing the spatial 
correlation (designated SC in the program) calculation. However, the program 
is designed strictly for panel-type response and only determines response per- 
pendicular to the plane of the panel. 

B.2 ORIENTATION OF PANEL 

In the SNAP analysis, the axis perpendicular to the plane of the panel 
must be of the first coordinate axis. The plane of symmetry, if used, must 
be the plane of the first and third coordinate axes. 
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B.3 MAIN PROGRAM LISTING 


' FOR , L I PSP,RSP 


1 

2 

3 

4 

5 

6 

1 

3 

3 

4 


I NTFGFP R 

DIMENSION RM ( 6 * 600 ) , RP(3,1200) , RC ( 3600 ) » S(546.10) , X(546> 

, V ( 546 ) , I L ( 546 ) , A(546) * F(500) , DS(500) * AS<50C) 

« WS(10) , l M ( 10) , WSQ (ICO) , GMAS(IOO) , D(10) 

« P ( 500 ) , PR ( 500 , 2 ) * VS (2*2) » TL(6) , HD < 6 , 2 ) 

♦ HG 1(1?) , HG2 (12) , !P(2) , HG(12,2) ♦ TX(lO) 

« TY (10) , TLM (2*2) * BLM(2.2> , GLM ( 2 ) * RLM(2> 

« SR ( 1 0 ) 

EQUIVALENCE (RM,RP,RC,F) « ( RC ( 50 1 ) , DS , PR ) * ( RC ( 1 00 1 ) , AS ) 

, ( RC ( 200 1 ) » A ) ♦ (RC(3),JT> 

, ( RC ( 4 ) »MTi , <RC(5),WSG) , ( RC ( 1 05 ) , GMAS ) 

* ( RC ( 1 50 1 ) , P ) * ( VS , VD ) , (VS(1«2)*VA) 
, (HG,HG1) , ( HG ( 1 ♦ 2 ) * HG2 ) 

NAMELIST /JNTS/NJ, IL*NM* IM* JTNiRtNS 


1 /FREQ/NF »F«AtD*P( TL 

DATA P I /3 . 1 4 1 59265/P2 I /6 . 283 1 B53/P4 1/12. 5663706/ I P/0 , 1 /CBF/ *00000 
lO*/ HD/72HD I SPLACEMENT RESPONSE IN INCHES ACCELERATION RESPO 

2NSF IN G*S / 

3 HG 1/72HTPS PANEL VA 

4R = P. E / 


JU 1 N I 


5 HG2/72H 

6S = O « E 

7TY/ * 

ATX/' 

READ (5.JNTS) 

WRITE (6.JNTS) 

IF (NS.LT.O) 

IF (NS.GT.NJ) 

N7 = NS + 1 
CALL NTRAN (10,10) 

CALL NTRAN ( 10, 2,204, RC, LSTA) 

CALL NTRAN ( 10,22 ) 

IF (L C TA ) QOOO 

WRITE (6,2) RC( 1 ) , RC ( 2 ) ,JT,MT, ( WSQ (I),I=1,MT),( GMAS (I)»I=1»MT) 

2 FORMAT ( • 1 DATE *,A6,» , TIME *,A6»* * JT = *,I6«* , MT = *,16/ 


/B2M/* 00 */BLM( 1,1) ,BLM( 1 ,2) , GLM ( 1 )/3*0.0/ 

POWER SPECTRAL DENSITY */ 

FREQUENCY (HZ) •/ 


NS = O 
NS = NJ 


1 ( 1 0X,5E20.8 ) ) 

IF (MT-IM(NM)) 9P99 
IF (JT-JTN) 9999,, 9999 
J3T = 3* JT 
J6T = 6* JT 
DO 4 1=1, NM 

J = !M( I ) 

WS ( I ) = WSQ ( J ) 

4 CONTINUE 

CALL NTRAN ( 1 O , 2 , J3T , RC , LSTA ) 
CALL NTRAN (10,22) 

IF (LSTA) 9000 
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no ft K= 1 ,NJ 
J » ! I (O 


X ( (O = PP ( 2 » J ) 

Y(K> = PP (3* J) 
ft CONTINUE 

CALL NTRAN (10«7,IM(1>) 

DO IP M=1«NM 

CALL NTRAN ( 10,2, J6T , PC »LSTA ) 
CALL NTRAN (10,22) 

IF (LSTA) 9000 
DO 10 J=1,NJ 
1 = I L ( J ) 

S( J,M) = RM (1,1) 


10 CONTINUE 

SR ( M ) = RM ( 1 , R ) 


IF (NM-M) 12*12 

I = IM(M+1) - IM(M) - 1 

IF ( I ) 12,12 

CALL NTRAN ( 10,7, I ) 

12 CONTINUE 

READ (F,FRFO) 

DO lft M=1 ,NM 

D(M) = 4.0*D(M)*D(M) 
lft CONTINUE 

DO ftO L= 1 , NE 


DO 

DO 


DO 


DA = 0.0 
W = P2 I *F ( L ) 

W 2 = W*W 
WDC = W/ 13400*0 
FO M= 1 ,NM 
AJ = 0.0 
IF (NS) ,41 





IF 


40 J=1 , NS 

AJ = AJ + (S( J,M)*A( J) )**2 

NK = J - 1 
(NO ,40 
3D K= 1 , NK 

DL = SORT < (X( J)-X(K) )**2 + { Y ( J ) - Y ( O ) **2 ) 

SC = SIN (DL) / DL 

AJ = AJ + SC*S( J,M)*S(K,M)*A ( J)*A(O*2.0 


* WDC 



30 CONTINUF 


40 CONTINUE 

IF (NS-NJ) ,4ft 

41 CONTINUF 

DO 4S J=N7,NJ 

DL = 2.0 * X(J) *WDC 
SC = SIN (DL) / DL 

AJ = AJ + 0 • *5 * (1.0+SC) * ( S( J,M)*A( J) )**2 
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4? 
4 0 


(NS) 

.4 

0 





4? 

K = 

1 . NS 






OL 

= 

SORT 

( ( X ( J ) 

-X ( 

K) 

) **2 

+ (Y( J)-V(K) )**2) 

SC 

= 

SIN 

(DL ) / 

OL 




AJ 

- 

AJ + 

SC*S ( J 

.M) 

*S 

(X.M) 

*A( J)*A(K)*2.0 


CONTINUF 

CONTINUF 

NX = J - 1 
IF (NX-N7) 45 

DO 44 K=N7.NX 


* WDC 


DL = 
SC1 = 
DL = 
SC? = 
SC = 
AJ = 


SORT 

( (X( J)-X(K) ) **2 + 

( Y ( 

J)-Y(K) >**2) 

* WDC 

SIN 

(DL) / DL 




SORT 

( (X( J)+X(X) )**? + 

( Y( 

J)-Y(K) )**2) 

* WDC 

SIN 

( DL ) / DL 





SCI + SC? 

AJ + SC*S( JiM)*S<KtM)*A( J)*A(X) 


44 CONTINUF 

45 CONTINUE 
45 CONTINUE 

DA = DA + AJ* (SR(M)**2) / ( ( WS ( M ) “W2 ) **2 + D < M ) *WS ( M ) *W2 ) 
50 CONTINUE 

OS(L) = P(L) * DA 

AS ( L ) = 0.6709E-5 * W2 * W? * DS(L> 

60 CONTINUF 

FI = F(?) - F( 1 ) 

FIV = F ( NF ) - FINF-l ) 

VD = DS( 1 )*FI + DS(NF)*FIV 
VA = AS ( 1 ) *F I + AS ( NF ) *F I V 


64 

L = 0 . NF 



FI 

= F ( L ) 

- F(L- 

?) 

VD 

= VD + 

DS ( L— 1 

) *F I 

VA 

= VA + 

AS (L-l 

) *F I 


NOT 


reproducible 


64 CONTINUF 

VD = O.S*VD 
VA = 0 . 5*VA 
VS( 2.1 > = SORT (VO) 

VS (2,?) = SOPT(VA) 

'a/RITF (6.6 C) TL.R.VS. (DS( I ) « 1 = 1 » NF ) 
66 FORMAT (»1TPS PANEL *.6A6»6X.*J0 
1 E * « 6X « • V AR = • .E16.8.6X. *RMS = ».E1 


2SP0NSE' »6X. • VAR = * * E 1 6 • 8 . 6X . » RMS = 
OFPECTPU^ • / ( 'O' 5F?0 . 0 ) ) 

WRITE (6.6 8) (AS(I).Isl.NF) 

68 FORMAT ( • 1 ACCELERAT I ON SPECTRUM* /( 
WRITE (6.70) ( F( I ) . I = 1 .NF) 

70 FORMAT ( • lFREQUFNCIES»/( ' 0*5E20»8) 
DO 7 p 1=1.6 

HG1 ( 1+2) = TL ( I ) 


INT * . 15/* ODISPLACEMENT RESPONS 
6.8.* inches*/*oacceleration RE 

'.E16.8.4H G»S/' ODISPLACEMENT 


• 0'5E20.8) ) 
) 
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72 CONTINUE 

CALL PLOTIN (9) 

CALL BNBCDV( ( FLOAT ( R ) +0 . 0 1 ) ,NPR,NDS) 

NPR = OR ( NPR « CRF ) 

K 1 = 46 - NDS 

CALL MOVST ( NDS «K1 i HG2 iO ( NPR ) 

RLM (1) = E ( NF ) 

CALL MXMNLG (NF.DS,BLM( 2, 1 ) ,TLM( 1 , 1 ) ) 

CALL MXMNLG ( NF « AS « BLM ( 2 , 2 > , TLM ( 1 » 2 > > 

CALL LOGLIM < 1 . 0 , F ( 1 ) , F ( NF ) * GLM ( 2 ) , RLM C 2 ) ) 

CALL LOGLIM ( 1 • 0 « BLM ( 2 « 1 ) * TLM (1,1), BLM (2*1 ) « TLM ( 2 * 1 ) > 

CALL LOGLIM (1,0, BLM (2*2) » TLM (1*2), BLM (2*2) » TLM (2*2) ) 

no on k=i ,2 

OO 78 1=1,6 

HG2( I ) = HD ( I ,K ) 

7P CONTINUE 

OO 80 1=1,? 

CALL BNBCDV ( VS {!»(<)♦ NPR « NDS > 

NPR = OR ( NPP « CBF ) 

CALL MOVST ( B , 60 , HG ( 1 « I) ♦ 0 ,NPR ) 

CALL BNBCDV ( ( FLOAT ( NDS ) ) ♦NPR.Kl ) 

NPR = OR (NPR, CBF) 

CALL MOVST C 3 « 67 , HG ( 1 ♦ I > , 1 * BZM ) 

IF (NDS) ,80 
K? = 70 — K 1 

CALL MOVST (K1 ,K?,HG( 1 , I ) , O , NPR ) 

IF (NDS) ,,80 

CALL MOVST < 1 , 67 , HG ( 1 , I ) , 5 , BZM ) 

80 CONTINUF 

DO 82 1=1,2 

CALL GR I DPL ( I P ( I ) , I P ( I ) , 20 , 1 , NF , F , PR ( 1 , K ) , GLM ( I ) , RLM ( I ) , BLM ( I , K ) , 
1 TLM ( I , K ) , TX , TY , HG 1 , HG2 ) 

82 CONTINUE 
90 CONTINUE 

CALL ENDPLT 
GO TO 9999 

9000 WRITF (6,9001) LSTA 

QOOl FORMAT ('ONTRAN RFAD ERROR , LSTA = *,I4> 

CALL FRROFF 
OR99 STOP 
END 


not reproducible f! 
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B.4 SUBROUTINE PROGRAM LISTINGS 


•FnR,L! PTN , PTN 

SUBROUTINE PLOT IN (ID 

ry tern a i tamliv 


nIMFMMON flnA°Y ( 22 ) 

DATA ADARY/6HONF HA , 6HRDCOPY ♦ 20*6H / 

DATA DL/20.0/ 

CALL IDFNT (IC.ADARY) 

CALL SFTMIV (24,0,30,54) 

CALL CHFI7V (2,2) 

CALL RITSTV ( 1 2 , 1 8 , TABL 1 V ) 


2 


4 


8 
1 O 


rfturn 

ENTRY GRIDPL (LX»LY» IP* IL 1 N 1 X 1 Yt XL , XR » YB » YT » TX *TYiHl » H2 ) 
DIMENSION H 1 ( 1 2 ) , HP (12) ♦ TX(lO) , TY(lO) , X(N) , Y ( N ) 

CALL SMXYV (LX»LY ) 

CALL DXDYV ( 1 » XL , XR , DX , LVR . LVL , NX » DL , I ER > 

CALL DXDYV (2i YB » YT » DY , LHR « LHL ,NY , DL « I ER ) 

IF (LX) * 2 
DX = 1.0 
LVR = 10.0 
LVL = 10.0 
IE (LY) « 4 
DY = 1 .0 

LHR = 10.0 

LHL = 

CALL 
CALL 


10.0 
GRID IV 
R ITF2V 
CALL RITE2V 
CALL RITE2V 
CALL RITE2V 
ENTRY ADLINE (IP,IL,N,X»Y) 

CALL X55CLV 1 ( X ( 1 ) * NX 1,1 ER ) 

CALL YSCLV1 ( Y < 1 ) ,NY 1 , I ER ) 
CALL PLOT V ( NX 1 » NY 1 « IPiO ) 
DO 10 K = 2 , N 


( 4 , XL * XR * YB » YT , DX * DY * LVR , LHR , -LVL , 
( 148. 1 OIF. 1023.90, 1 , 72 , 1 , H 1 , I ER ) 
(148,991,1 023 , 90 , 1 . 72 , 1 , H2 . I ER ) 
(220,8, 1023,90, 1 , 60 , 1 , TX , I ER ) 
(8,140. 1023, 180, 1 ,60, 1 ,TY» IER) 


•LHL, NX, -3) 


CALL XSCLV1 (X(K) ,NX2, IER) 
CALL YSCLV1 ( Y ( K ) . NY2 , I ER > 
CALL PLOTV (NX2.NY2, IP.O ) 

IF (IL) 8,8 

CALL LINEV ( NX 1 , NY 1 , NX2 , NY2 ) 
NX 1 = NX 2 
NY 1 = NY2 
CONTINUE 
RETURN 
ENTRY ENDPLT 

CALL ENDJOB 
RFTURN 
END 












& 
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•FOR , L I MML » MML 

SUBROUT INF MXMNLG <N,X.B,T> 

DIMENSION X(N) 

t = x< n 

B = T 

DO 1 K=2«N 

IF (X(K) .LT. 1 .0E-38) X(<) = 1.0E-3B 

IF (X(K).GT.T) T = X ( K ) 

IF (X(K).LT.B) R = X ( K ) 

1 CONTINUE 
RETURN 
END 


•FOR, LI LGM,LGM 

SUBROUTINE LOGLIM ( F «B» T » BL , TL ) 

BL = ALOGIO <B> / F 

IF (BL.LT.O.O) BL = BL - 1.0 
BL = F*A I NT'( BL ) 

TL = ALOGIO (T) / F 

IF (TL.GT.O.O) TL = TL + 1.0 
TL = F*A I NT ( TL ) 

IF ( (TL-BL) .GT.7.0) BL = TL - 7.0 
BL = 10.0 ** BL 
TL = 10,0 ** TL 
RETURN 
END 
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B.5 DATA CARD INPUT 


This input data is read in Namelist format. Namelist JNTS contains 


JTN 

R 

NJ 

IL 


NS 


NM 

IM 


total number of joints in the finite element model 

joint number where response prediction is desired 

number of joints to be used in the double area integration 

array of NJ joint numbers to be used in the double area integra- 
tion (If the plane -of- symmetry method was used in SNAP, any 
joint numbers located in the plane of symmetry must precede 
the other joint numbers. ) 

number of joints in the IL array that are located in the plane 
of symmetry, if the plane -of- symmetry method was used in 
SNAP (If the plane -of- symmetry method was not used, NS 
must be equal to NJ. ) 

number of panel modes to be used in the analysis 
array of NM mode numbers to be used in the analysis. 


Namelist FREQ contains 


TL 

A 

D 

NF 

F 


36 character title that will appear on the plots and printout 

array of NJ areas that contain the panel area represented by 
the corresponding joint number in the IL array 

array of NM damping values for the corresponding modes in 
the IM array 

number of frequencies for power spectral density calculations 

array of NF frequencies where power spectral densities are 
to be calculated 


P 


array of NF values of the power spectral density of the acoustic 
excitation field for the corresponding frequencies in the F array. 
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B.6 SAMPLE DATA CARDS 


t.JNTF JTM 

= 50/4 

, 

R = 53 

, 

N8 = 

8, 


NJ = 96, 
IL = 5. 

7, 

9, 

1 1 , 

12, 

1 4 , 

16, 

18, 

47, 

49, 

51 , 

53 , 

54 , 

56 , 

58, 

60 , 

80, 

91 , 

93, 

95, 

96, 

98, 

100, 

102, 

131 , 

1 33, 

135, 

1 37, 

138, 

140, 

142, 

144 , 

173, 

1 75, 

177, 

179, 

180, 

1 82, 

184, 

186, 

215, 

217, 

219, 

221 , 

222, 

224, 

226, 

228, 

257, 

259, 

261 , 

263, 

264 , 

266, 

268, 

270, 

299, 

301 , 

303, 

305, 

306, 

308, 

310, 

312, 

34 1 , 

34 3, 

345, 

347, 

348, 

350, 

352 , 

354 , 

383, 

385, 

387, 

389, 

390, 

392, 

294, 

395, 

425, 

427, 

429, 

431 , 

432, 

434 , 

436, 

438, 

446, 

448, 

450 , 

452 , 

453, 

455, 

457, 

459, 

NM s 5 , 
5FND 

I M = 

1 , 2 

,3,4 

, 5 , 





5FRFO TL = 

36H 

HAYNES 25 ALLOY L 

A = 3*1.5, 

2*1.125, 

3*1.5, 

3*3.0 , 

2*2.25, 

3*3.0, 

3*3.0 , 

2*2.25, 

3*3.0, 

3*3.0, 

2*2.25, 

3*3.0, 

3*3.0 , 

2*2.25. 

3*3.0, 

3*3.0, 

2*2.25, 

3*3.0, 

3*3.0 , 

2*2.25, 

3*3.0, 

7*3.0, 

2*2.25, 

3*3.0, 

3*3.0 , 

2*2. 2F, 

3*3.0, 

3*3.0 , 

2*2. 2*, 

3*3.0, 

3*2.5 , 

2*1 .875, 

3*2.5, 

3*2.0 , 

2*1.5, 

3*2.0 , 

0 = 5*0.0 5 





NF = FI, 


F = 10, 20, 30, 

40, 50, 

60, 70, 

* 

o 

O' 

» 

c 

CO 

100, 106, 

110, 

1 1 3.09774 , 






116, 120, 

130, 140, 

150, 160, 

170, 180, 

185, 190, 

193, 

195, 197.7^882, 






200 , 

202, 204, 

206, 208, 

210, 211, 

212, 


212.65703, 213 

, 214, 

215, 




215.82162, 215 

.00176, 2 

15.O81Q0, 




217, 

218, 220, 

222, 224, 

227, 230, 

235, 240, 

250, 

P = 20.5F-5, 

18.2E-5, 

15.0E-5, 

1 2.9E-5, 

12.9E-5, 

16. OE-5, 

1 8.5F-=: , 5 

*20 • 5F— 5 , 





3*20. 5E-5, 



20 • OE— 5 , 

19.0E-5, 

17. 5E-5 , 

1 6.3E-5, 

14.8E-5, 

1 3.3E-5, 

1 2.5F.-5, 

1 1 . 8E— 5 , 

1 1 . 4E— 5 , 

1 1 . 1 E — 5 , 

10.7E-5, 

1 0.3E-5, 

10.1 F-5, 

9.9E-5, 

9.8E-5, 

9.6E-5, 

9.5E-5, 

Q.4F-5, 

2*9. 3E-5, 

9.2E-5, 

9. IE-5, 

9. 1 E — 5 , 

3*9 . OE— 5 , 

8.9F-5, 

8.8E-5, 

8.7E-5, 

8.6E-5, 

8.4E-5, 

8.2E-5, 

8.0E-5, 

7.6F-5, 

7.2E-5, 

6.5E-5, 


5FNO 
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